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Performance of a measurement-driven ‘adiabatic-like’ quantum 3-SAT solver 

Simon C. Benjamin 

Department of Materials, University of Oxford, Parks Road, Oxford 0X1 3PH, United Kingdom 

I describe one quantum approach to solving 3-satisfiability (3-SAT), the well known problem in 
computer science. The approach is based on repeatedly measuring the truth value of the clauses 
forming the 3-SAT proposition using a non-orthogonal basis. If the basis slowly evolves then there is 
a strong analogy to adiabatic quantum computing, although the approach is entirely circuit-based. 

To solve a 3-SAT problem of n variables requires a quantum register of n qubits, or more precisely 
rebits i.e. qubits whose phase need only be real. For cases of up to n = 26 qubits numerical 
simulations indicate that the algorithm runs fast, outperforming Grover’s algorithm and having a 
scaling with n that is superior to the best reported classical algorithms. There are indications that 
the approach has an inherent robustness versus imperfections in the elementary operations. 


I. SCOPE AND NATURE OF THIS DOCUMENT 

This document contains my notes on a series of nu¬ 
merical experiments I have performed to assess a quan¬ 
tum approach to solving classical 3-SAT problems. I use 
measurements to drive a register of qubits toward a target 
state satisfying a condition related to the 3-SAT (actually 
an instance of Bravyi’s quantum 3-SAT [1]). The target 
state can be slowly evolved from a trivial initial target 
to a final state that reveals the 3-SAT solution; this is 
analogous to adiabatic quantum computing [2] but using 
only circuit-like operations. 

The numerical results I present are unverified - they 
have not been independently reproduced, therefore either 
outright code errors or more subtle numerical inaccura¬ 
cies may be occurring. I have not previously studied 
3-SAT problems, therefore the instances that I generate 
may be the ‘wrong’ choice and may include too many 
cases that are easy classically - moreover I do not com¬ 
pare the performance of the quantum approach with a 
classical solver, I merely cite the reported performance 
bounds for such solvers. 

The purpose of this document is to serve as a basis for 
discussion. I try to introduce all the terms/ideas I use 
without presuming any specific expertise from the reader, 
beyond undergraduate quantum physics and a basic un¬ 
derstanding of algorithms. For experienced quantum re¬ 
searchers, the short version is given in the box. 


II. SATISFIABILITY 

Boolean satisfiability problems involve a ‘proposition’ 
P{xi,X2, ■■■Xn) = TRUE where the n boolean variables Xi 
may each take the value TRUE or FALSE. Typically the 
challenge is to determine whether any choice of values for 
the variables Xi can satisfy the formula. 

In the particular case called 3-satisfiability (3-SAT) 
the expression P{xi,X 2 , ■■■Xn) is formed by AND’ing 
together a number of clauses, and each clause involves 
three variables OR’ed together. 

Example: (xiV-'a; 2 Va: 3 )A(-'a:iVa:: 4 Va; 5 )A... = TRUE. 


Here the symbol A denotes the logical AND, V denotes 
the logical OR, while ^ denotes the logical NOT which 
is also referred to as negation. Table I summarises these 
standard elements. Note that verifying a candidate solu¬ 
tion merely involves checking each independent clause is 
TRUE; if-and-only-if they all are, the 3-SAT is satisfied. 


Overview of the paper 

I looked at an algorithm (actually a few variants of 
an algorithm) that seeks to solve 3-SAT problems 
by applying a series of post-selective measurements 
to a register of qubits. Each boolean of the 3-SAT 
corresponds to a unique qubit, and each quantum 
measurement involves checking that a clause of the 
3-SAT is satisfied. The extra trick is that the single¬ 
qubit states representing TRUE and FALSE are not 
in general orthogonal. 

The process of performing a general clause check is 
shown in the box on page 4. There is a unique state 
that must pass all the clause checks (Eqn. 2) if the 
3-SAT has a unique solution (and a corresponding 
subspace if multiple solutions, see Appendix). 

I’ve looked at a variant where the basis of 
TRUE/FALSE slowly evolves until becoming or¬ 
thogonal, this is somewhat analogous to adiabatic 
quantum computing - see the flow diagram Fig. 1 
and data in Figs. 2 & 3. I’ve also looked at a 
rather simpler variant where one creates (‘sculpts’) 
the satisfying state for a fixed basis, before measur¬ 
ing it for partial information on the 3-SAT solution 
- see flow diagram Fig. 4 and Figs. 5 & 6. This 
latter approach is the one I’ve studied more fully 
using numerical experiments. I’ve also looked at a 
hybrid, see Figs. 7 and 8. 

Some of these approaches outperform both Grover’s 
and leading conventional algorithms, over the very 
limited range I can simulate and with the caveat 
that I can’t specify a fixed time within which the 
solution will certainly be found. 

I have not properly investigated error tolerance, but 
I have tried adding 1% or 2% random over/under 
rotation to the simulation, and this seems well tol¬ 
erated, see Fig. 10. 
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TABLE I. Table defining the AND, OR and NOT. 

III. 3-SAT PROBLEMS CONSIDERED HERE 

In the analysis described here, I generate 3-SAT in¬ 
stances randomly but in such a way that 

1 . All clauses are distinct from one another. 

2. Every variable Xi occurs at least once in the 3-SAT 
proposition in positive form Xi, and at least once 
in negated form -^Xi. 

3. All clauses involve three different variables, i.e. nei¬ 
ther (xi V xi V X 2 ) nor (xi V -^xi V X 2 ) would be 
generated. 

4. Given n variables I generate = round(i? n) 
clauses, with R = 4.267. 

Having generated a 3-SAT formula in this way I then 
determine the number of solutions. Ns- This is a trivial 
task classically for the problem sizes I consider. Gener¬ 
ally I will be wish to have 3-SAT problems with a specific 
Ns for testing the quantum algorithm; if the randomly 
generated 3-SAT has some other Ns I discard it and gen¬ 
erate a fresh formula. For the majority of data I present 
here I have chosen Ns = 1, i.e. problems with a ‘unique 
satisfying assignment’ (USA), however I do also consider 
cases with Ns > 1. While it is quite rare to have exactly 
one solution, so that these 3-SAT instances are special, 
I choose this to reduce variably in the challenge of solv¬ 
ing them. Note that criterion ( 1 ) is simply to ensure 
proper counting of the number of clauses - if two clauses 
were identical, one could be dropped without changing 
the 3-SAT problem. Similarly (2) is imposed because 
otherwise the effective number of variables is reduced: 
for example if a variable only occurs in positive form 
Xi without —'Xi appearing anywhere in the formula then 
trivially one chooses Xi=TRUE and all clauses involving 
Xi become TRUE. Meanwhile (3) is merely for simplic¬ 
ity - in principle recurrence of a variable within a clause 
is legitimate for 3-SAT but I do not consider such cases 
(although there is no apparent difficulty with such a gen¬ 
eralisation). Point (4) is to maximise the chances that a 
randomly generated set of clauses will have a small but 
non-zero number of solutions Ns- I rely on the results 
of previous SAT studies which have found that there is 
a key threshold in the ratio R between the number of 
clauses and the number of variables - instances below 
the threshold are likely to be satisfied with multiple so¬ 
lutions, and instances above it are unlikely to be satisfi- 
able. This threshold is still somewhat debated, and only 


becomes sharp as n —> 00 , but the value 4.267 has been 
suggested [3]. 

IV. CLASSICAL 3-SAT SOLVERS 

The task of efficiently solving 3-SAT instances is a long 
studied challenge in computer science. It is believed that 
any algorithm will require a running time which is expo¬ 
nential in n [4]. Clearly the problem can be solved by a 
brute force search over all 2" variable assignments. Then 
each failed test of a potential assignment prior to the so¬ 
lution will require some portion of the W clauses to be 
checked, so that t = ^ 2 ” is a trivial upper bound to 
the number of clause checks required. Throughout this 
paper I take the expected number of clause checks 
to be the metric for the running time of an algorithm; 
this is fair only if the time required for a clause check is 
independent of the size of the problem, of course. 

The naive classical bound can be dramatically im¬ 
proved upon. A series of publications have succeeded 
in reducing the base K in the the expression for the ex¬ 
pected running time, t oc if". The current record holder 
appears to be an algorithm due to Paturi et al [5] which 
boasts K = ^ 1 ^ 307 . 


V. A QUANTUM 3-SAT SOLVER 

A simple application of Grover’s search algorithm 
would yield a quadratic improvement over the brute force 
classical search, from t oc 2" to t oc (v^)" i-e. K = 1.41, 
but this is in fact inferior to the best classical solutions 
as noted above. It has been suggested [ 6 ] that one of 
the high performing classical algorithms can be acceler¬ 
ated by replacing the random search with a Grover-type 
coherent evolution, leading to a hybrid algorithm with 
running time oc (1.153)". 

Here I describe a different quantum approach based on 
repeatedly making projective measurements, each corre¬ 
sponding to evaluating the truth value of a generalised 
clause. I begin by assigning a qubit to represent each 
boolean variable Xi in the problem. The general state 
of an isolated pure qubit is of course |4>) = a|0) -I- /3|1), 
where a and /3 can be complex, but in this analysis I 
find it suffices to consider qubits on the real great circle 
of the Bloch sphere, i.e. where a and /3 are real num¬ 
bers. Qubits of this kind are sometimes called rebits. In 
the following I choose my basis so that, if I am using 
an orthogonal basis for representing truth values, then 
the state |1) will be identified with boolean TRUE, and 
|0) will be identified with FALSE. (Obviously, rotating all 
states and bases in what follows will yield the same algo¬ 
rithm with, say, |-|-) and |—) as the truth values, if that 
is preferable). In general I consider cases where TRUE 
and FALSE are not orthogonal. 

I begin by describing a naive all-orthogonal approach 
to establish the framework. Consider a operation called 
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“Q-NOR^”, which involves 4 qubits, a, b, c, d. The 
operation simply applies a flip to qubit d if, and only if, a, 
b and c are all in the |0) state; it does nothing otherwise. 
By ‘ flip’ I refer to the Pauli x-operator. Therefore this 
four-qubit operation is merely a control-control-control- 
not gate acting when the controls are zeros. 


TABLE II. Q-NORz operation 


a 

b 

c 

action on d 

|0) 

|0) 

|0) 

flip: |0) ++ |1) 

|0) 

|0) 

ll) 

none 

|0) 

ll) 

|0) 

none 

ll) 

ll) 

ll) 

none 


the succession of clause checks, and so the probability of 
success is simply the amplitude-squared on the correct 
term, 2“"’. Therefore if we seek to solve our 3-SAT in 
this way, we expect to have to restart ~ 2” times before 
finding a solution, equivalent to the worst case of the 
brute force classical approach. 

Now I describe a simple generalisation of this setup: a 
quantum algorithm that is again based on repeated mea¬ 
surement of the truth value of clauses, but within a basis 
structure such that TRUE and EALSE are not generally 
orthogonal. As a segue to that general case, first suppose 
that the basic multi-qubit operation that our computer 
can directly perform is not in fact “Q-NOR^” but rather 
“Q-NORx” as shown in the following table, which uses 

|-)^(|0)-|l))/v^. 


If I prepare d in state |1) and apply the Q-NORz op¬ 
eration, and then measure d in the z-basis (i.e. the basis 
with outcomes |0) and |1)), then the outcome tells me 
the truth value of the a V 6 V c clause: if-and-only-if a, 
b and c are all FALSE then d will be measured as |0), 
denoting FALSE. If I wanted to test the clause aV ^bV c 
then I would simply invert b before this process (i.e. ap¬ 
ply a x—gate to 5), and again afterwards. What is the 
effect of this kind of process when the qubits a, b and c 
are in some general state? If we fail the clause check, by 
measuring |0) on the ancilla, then we project a, b and c 
onto state |000) (where I now write all three qubit states 
within a single ket for compactness). Conversely if we 
pass the check by measuring the ancilla in state |I), then 
we have projected out (removed) any component |000) 
that had been present in the state of a, b, c. 

Now consider the effect of a series such clause checks, 
one for each clause in our 3-SAT problem. Note that 
all such clause checks will commute with one another, 
because they are simply projections in the z-basis. Con¬ 
sequently if all clause checks are successfully passed then 
subsequently measuring the qubits in the z-basis will pro¬ 
duce a satisfying assignment of truth values for the 3-SAT 
proposition. 

Suppose that we initialise our qubits such that each 
is in the ‘plus state’ |-|-) = (|0) + |I))/-\/2, i.e. an equal 
superposition of our orthogonal TRUE and FALSE states. 
Then we perform each clause check, one after another, 
until either we see a ‘failed’ outcome on the ancilla, or 
else we have passed one clause check for each clause in 
our 3-SAT problem. If indeed we succeed then we can 
directly measure the measure the qubits one at a time 
to obtain a satisfying solution; what is the probability 
of such a success? To answer this note that our initial 
state is simply the equal superposition of all states in the 
z-basis: 


2-5(|0...00) + |0 ...01)+ ... +11 ...11)). 

If we are dealing with a 3-SAT that has a single solution, 
then all but one of these terms will be projected out by 


TABLE III. Q-NORx operation 


a 

b 

C 

action on d 

I-) 

H 

|-> 

flip: |0) ++ |1) 

I-) 

I-) 

1+) 

none 

I-) 

1+) 

I-) 

none 

1+) 

1+) 

1+) 

none 


Well we can simply recover the former clause check 
process by individually rotating the three control qubits 
appropriately, both before and after the application of 
“Q-NORx”. We can use a T(|) gate, where 


Y{9) = 


cos I sin I 
— sin I cos I 


so Y { - ] = ^ 


(f) 


^2 I -1 


The T(|) rotation maps state |0) to |—) and maps state 
|1) to |-|-). Therefore applying this operation, separately 
and independently, to the three qubits a, b and c will 
map from any of the input rows in Table (II) to the corre¬ 
sponding row in Table (III). As before we prepare ancilla 
d in state |1) prior to applying the Q-NORx operation, 
then subsequently measure d in the z—basis and say that 
the clause check is ‘passed’ if-and-only-if the measured 
state is |1). Finally we apply the inverse rotation Y(—^) 
to each of the qubits a, b, and c; this returns them to 
their original basis. 

Recall that in the naive protocol if we have a clause in 
which one, or more, of the variables appears in negated 
form then we should invert the corresponding qubit (s) 
|0) ++ |1) before, and then after, the operation. But 
now that we are in any case applying rotations to our 
qubits, we notice that we can simply apply Y{—^) prior 
to the Q-NAND^ (rather than T(|)) to any qubit whose 
corresponding boolean appears in negated form in the 
clause. This maps |0) to |-|-) and correspondingly |I) to 
(—1)|—), where we note an extra phase of —I in the latter 
case. After the ancilla measurement we of course apply 
Y{^) to reverse the former rotation, and reabsorb the 
phase. 
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Thus our earlier protocol for clause evaluation is 
equivalent to the following steps, if we choose 0 = |. 
I refer to this process as a clause check operation. 


Clause check operation: 

1. Select a clause; let the three boolean variables 
to which it refers be Xi, Xj, Xk, and denote the 
corresponding qubits as Qi, Qj and Qk- 

2. Rotate qubit Qi by Y(0) if variable Xi appears in 
positive form, or by Y{—9) if the variable appears 
in negated form. Rotate the other qubits by the 
analogous rule. 

3. Apply Q-NORx (Table III) targeting an ancilla 
prepared in | 1 ), and measure the ancilla in the | 0 ), 
| 1 ) basis. Abort if |0) is seen. 

4. Reverse all three rotations applied in Step 2. 


If indeed 0 = f then this is entirely equivalent to the 
naive protocol, with the same low success probability i.e. 
we have probability 2“" of ‘passing’ all the clause checks, 
when we start from an initial state | + + • • • +) and there 
is exactly one solution to the 3-SAT problem. 

Generally when we pass a clause check operation, this 
leaves the register in a state with zero projection on the 
particular three-qubit state that would have failed (i.e. 

the state that would have rotated to |-)). Asking for 

the state that passes all clause checks, for some particular 
value of 0, is therefore a special case of the quantum k- 
SAT problem introduced by Bravyi [I]. 

It is helpful to think of a clause check process as check¬ 
ing the truth values in a basis defined by 9. Suppose that 
a particular clause contains a boolean a for which the cor¬ 
responding qubit is Qa- If the variable appears in direct 
form, then classically the clause would certainly evalu¬ 
ate to TRUE if boolean a is itself TRUE. The equivalent 
statement for the quantum case is this: the clause check 
will certainly be passed if qubit Qa is in state y(—0)|-|-). 
The reason is simple: a clause check process begins by 
applying rotation Y{9) to Qa, thus leaving it in state |-|-), 
and then it is impossible that the ancilla will flip during 
the Q-NORx operation regardless of the other two qubits 
associated with the clause. Equivalently, if the clause con¬ 
tains the variable in its negated form then classically 
the clause will certainly be TRUE is a is FALSE, while 
the quantum clause check will be passed if qubit Qa is in 
state T(0)|-|-). The mapping is therefore, 

TRUE —>• y(—d)|-|-) = cos^|0) -I- sin 011) 0 = ^ ^ 

FALSER r(6»)|-h) = sin0|O)-bcos0|I). (I) 

We see that in the limit 0 = f we have the orthogonal 
states |1) and |0), whereas in the limit 0 = 0 the truth 
values are degenerate: both correspond to |-|-) (so that 
passing all clause checks is inevitable for an initial state 
|-|- -I- ..-f) but this is meaningless). 


For intermediate values of 9 the states are non- 
orthogonal: We cannot use a measurement to determine 
which state a qubit is in with certainty. Moreover the var¬ 
ious check operations need not commute', if two clause 
checks do not commute, then the order in which two 
clause checks are performed affects the resulting state of 
the qubit register. Two checks will not commute if they 
share at least one variable in common, and either (a) the 
two checks are using different 9 values, or (b) they are 
using the same 9 value but a variable appears in direct 
form in one clause and negated form in the other. That 
is, they do not commute if there is a conflict in the basis 
in which they are performing their projections. 

Consider the significance of non-commuting clause 
checks. If a clause check C is performed and is success¬ 
fully passed, and then that same clause check C is imme¬ 
diately repeated, then the second instance must succeed 
(since the first check projected out, i.e. removed, any 
part of the overall state that could fail). This remains 
true if one or more additional clause checks which com¬ 
mute with C are performed in the interval between the 
two instance of C. However if a non-commuting clause 
check V is performed between the two instances of C, 
then the second instance of C is not guaranteed to suc¬ 
ceed - the act of passing though clause check V will (in 
general) have re-introduced a finite probability of failing 
C. Now consider applying a full cycle of clause checks, i.e. 
one clause check for each clause of the 3-SAT, applied in 
some fixed order^. The fact that the system’s state has 
successfully passed one such complete cycle does not im¬ 
ply it would necessarily pass a second complete cycle - in 
fact it may fail the very first clause check of the second 
cycle, since several non-commuting checks will have been 
applied since that check was previously passed. 

Suppose that we do indeed apply a full cycle of clause 
checks. Given that the 3-SAT problem has exactly one 
solution, then there is one special state of the qubit regis¬ 
ter which is guaranteed to pass all checks. This is simply 
the following state, where I use the bold symbol in |0) to 
denote the entire n—qubit register: 

|0) = (g) y(L,0)|-f). (2) 

i 

where Li = —I if the satisfying logical value of the 
boolean variable Xi is TRUE, and Li = -|-1 if that value 
is Xi is EALSE. Any clause check will then rotate at least 
one of its three qubits to the |-|-) state, and therefore 

that check will succeed (since only the |-) row in 

Table III leads to a clause check failing). A converse ar¬ 
gument shows that only this state can satisfy all clause 
checks with certainty^, and this observation is verified 


^ In all cases I consider, I randomly select the order of the 3-SAT 
clauses at the start of an algorithm and then maintain that order 
through the algorithm 

^ To see the uniqueness of this state, consider the 2 " states of the 
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for every Ns = 1 problem for which I’ve performed a 
numerical simulation (hundreds of thousands of cases). 

This raises some interesting possibilities. If the state 
given in Eqn. 2 can be created, for any 9 > 0, then the 
qubit register encodes the satisfying solution to the 3- 
SAT problem. However in general one cannot determinis¬ 
tically obtain this information by directly measuring the 
register qubit-by-qubit. This is possible only in the limit 
6 — where one can simply measure in the |0), |1) basis 
to determine each boolean value FALSE or TRUE; but 
of course, simply setting 9 — initialising to |-|- -f ..-I-), 
and performing clause checks will lead to the low prob¬ 
ability of passing all clause checks. Alternatively 
we could start by performing clause checks with some 9 
close to zero (and which will therefore succeed with high 
probability) and then slowly increase 0, so as to gradu¬ 
ally transform the qubit register toward the desired final 
state with 0 = |. We might hope that the probability 
of successfully reaching that state is much greater than 
2 “”; this would be analogous to adiabatic quantum com¬ 
puting, where the parameters of the physical system are 
slowly changed and one hopes that the system will re¬ 
main in its ground state even as the nature of that state 
changes. 

An “adiabatic-like” approach 

A small change in the parameter 9 corresponds to a 
small change in the state \6) given by Eqn. 2, in fact 
|(6/'|0)|2 « 1 _ |g2 jf e = - 6» is small (c.f. the Zeno 

effect), and we may therefore hope to be able to guide 
the system to the ideal state \6 = ^) and arrive there 
with a probability that approaches unity for sufficiently 
slow evolution of 9. This intuition will prove to be cor¬ 
rect, but with the significant caveat that the probability 
of failure can only be made small for extremely slow evo¬ 
lution. We can select any rate of change of 9 simply by 
introducing an appropriate number of clause checks - in 
my simulations of this approach I set parameter 0 to a 
given value, perform one complete cycle of clause checks, 
and then increment 9 and repeat. See Fig. 1. 

Suppose that we try increasing 0 as a linear function 
of time, i.e. incrementing 0 by a small fixed amount <5 
after each cycle. What is the optimum value of SI One 
might imagine that it is best to make S as small as is 
necessary to make the probability of success approach 
unity. This intuition proves to be incorrect. To under¬ 
stand why, it is instructive to look at a simulation of a 
particular problem; in Fig. 2 I show data for a typical 


form of Eqn. (2), one for each of the possible sets of boolean 
assignments, i.e. each possible set of Li values. These states are 
independent, and therefore although they are not orthogonal for 
general 8 nevertheless they suffice as a basis to express any state 
of the qubit register. Decompose any general state in this way 
and then note that each element will fail at least one clause check, 
except for the element corresponding to the boolean assignment 
that satisfies our single-solution 3-SAT problem. 



FIG. 1. Flow diagram for the basic ‘adiabatic-like’ algorithm. 


USA 3-SAT problem involving 24 booleans. We see that, 
as expected, a slower evolution leads to a better chance 
of reaching the desired 0 = f final state. 

However, in fact one should not evolve more slowly 
than a certain optimum, and this optimum corresponds 
to a success probability that is still far below unity. This 
is because there is a peak to the failure probability when 
0 is already substantial, around in the mid-point in the 
evolution as can be seen from the middle graph of Fig. 2. 
The lowest graph in Fig. 2 shows that doubling the num¬ 
ber of increments from 92 (green) to 184 (blue) causes 
the success probability for a given ‘run’ to increase. How¬ 
ever the expected number of steps-until-failure on a each 
failed run increases by a similar factor because the peak 
probability of failure is now occurring after about twice 
as many cycles. The net consequence is revealed in Fig. 3 
where we see that going from green to blue corresponds 
to a slight increase in overall expected number of clause 
checks required before a 3-SAT solution is found. 

It is interesting to compare the performance of this 
simple measurement-driven algorithm to a direct use of 
Grover’s search algorithm. For the particular 3-SAT 
problem used for the simulations depicted in Figs. 2 and 
3, the overall performance with the optimum choice of 92 
full clause check cycles, is somewhat inferior to Grover’s 
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algorithm. The expected total number of clause tests in 
Grover’s would be 288,252 (see Appendix I), as com¬ 
pared to 386,462 for the present algorithm. If how¬ 
ever we alter the ramping function we can quickly im¬ 
prove the pe rforma nce of the present approach; by us¬ 
ing 9 = ^y'c/ctot, for the cycle, we reach an ex¬ 
pected number of clauses of only 190, 885 (this root-t 
evolution is shown later in Fig. 9). There are a number 




0 20 40 60 80 100 120 140 160 180 

Number of complete clause cycles 


FIG. 2. Data from a simulation of the ‘adiabatic-like’ ap¬ 
proach (Fig. 1) with a 3-SAT problem of 24 booleans and a 
unique solution. The red, green and blue traces correspond 
to three different rates of linear evolution of the parameter 
9. The slower the evolution, the higher the chance of reach¬ 
ing the desired state ultimate state \6 = ^); the fidelity with 
respect to this target is shown by the grey curves in the up¬ 
permost figure. However it is not optimal to go very slowly, 
as shown by the cost analysis in Fig. 3. 


of considerations that mean this comparison to Grover’s 
algorithm has limited meaning - for example, the per¬ 
formance of Grover’s is the same for all 3-SAT prob¬ 
lems of a given size and with a given number of solu¬ 
tions, whereas for the present approach there is variation 
(as discussed presently) and the data depicted in Figs. 2 
& 3 is merely typical rather than definitive. Neverthe¬ 
less is it interesting that the performance of Grover’s 
algorithm can be matched, roughly speaking, since the 
present measurement-driven approach is apparently very 
different. Moreover by adjusting the present approach its 
performance can be considerably improved. 

Although there is an analogy between the slow evolu¬ 
tion of 9 used here and the slow change of a Hamiltonian 
in AQC, already we are exploiting a feature which AQG 
does not offer: the moment that a run has failed, this is 
known and we can restart immediately without ‘throw¬ 
ing good money after bad’. Indeed, this is the reason 
that the square-root ramping algorithm outperforms the 



FIG. 3. A cost analysis for the same problem considered in 
Fig. 2. Curves show the expected total number of clause 
checks that are performed during runs that ultimately fail, 
assuming that we try enough times to achieve one success 
(and therefore reveal the solution to the 3-SAT). The curves 
are cumulative, showing the contribution to this total made 
at different stages (so, even though the cost of a failure late 
in the evolution is high, the contribution to the expected to¬ 
tal is small because such failures are rare). By adding this 
expected total to the number of clause checks in one success¬ 
ful run (96 x 102 = 9, 792 for the green trace), we can find 
the expected number of clause checks required to solve this 
3-SAT. This total of 386, 462 is optimal for a linear ramping 
of 6, but other ramps are better; see Table ?? and Fig. 9. 
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linear one - it ‘front loads’ the failure probability more 
aggressively, thus allowing us to fail faster! 

There are a number of approaches one might consider 
to refine this adiabatic-like approach, including more 
complex evolution of 6 over the course of the process 
or even attempting to ‘repair’ a state after a clause check 
fails. However I will not consider these further here, but 
instead I focus on simply reducing the expected num¬ 
ber of steps in a failed run by ‘front loading’ the failure 
probability even more aggressively. This approach will 
move away from the idea of an ‘adiabatic’ evolution, since 
we will now see a substantial failure probability for each 
clause check early in the process. 

A “sculpting” approach 

Consider an approach where we initialise all qubits in 
|-l-), select some fixed value of 9 (well above zero, yet well 



FIG. 4. The ‘sculpting’ algorithm aggressively ‘front loads’ 
the failure probability. 


short of the ideal 0 = f state) and simply begin perform¬ 
ing clause checks. The clause checks early in the process 
will be the most likely to fail, and (if passed) will have 
the greatest effect on the system’s state - somewhat as 
a sculptor might begin by carving away large pieces of 
stone. Our aim is to create state |0) to a high fidelity 
after a sufficient number of “passed” clause checks. Be¬ 
cause of the non-orthogonality of the clause check oper¬ 
ations, it may require many complete clause check cycles 
to reach high fidelity. As I presently discuss, although 
this target state is not the ideal 0 = f state, we can 
extract useful information from it and then repeat, as 
shown in Fig. 4. Alternatively, having achieved the state 
|0) we could then apply an adiabatic-like evolution as a 
kind of secondary phase, in which we guide the system 
onward from this 9 to the ideal as shown in Fig. ??. 


Since there is a finite projection between the initial 
state and the target state \6), we expect to to arrive at 
this target with at least probability 

p= |(6| + -I-1-)|^ = ^cos 

However since we are selecting a 9 value that is substan¬ 
tially greater than zero, the process is not an adiabatic 
evolution and we should not expect to achieve a success 
probability that is significantly higher than this bound. 
The simulations I have performed verify this. 

It is interesting to monitor the fidelity of the qubit reg¬ 
ister with respect to the target state |0) as the algorithm 
proceeds, i.e. as each successive clause check operation is 
successfully passed. Specifically, it is useful to note how 
many clause checks must be passed before the fidelity 
reaches a given threshold (since the fidelity only asymp¬ 
totically approaches unity, we must select some threshold 
below unity as a target). Having chosen a target fidelity, 
one can ask whether different 3-SAT problems take longer 
to reach the target and moreover how this varies as the 
number of booleans increases. 

Figure 5 shows AhiPid, the total number of successful 
clause checks that must be performed when we start from 
the initial state |-|- -I- ..-!-) and aim to reach |0) with a 
fidelity on the target state of 0.999. This is, A^hiPid is 
number of clause checks until 



|(e|4')|2 > 0.999 

where Idt) is the instantaneous state of the qubit regis¬ 
ter at some point during the algorithm. The interesting 
observation is that there is a spread in ApiPid even for a 
fixed choice of n and 9. Note that all problems generated 
for these simulations had a unique satisfying solution. 

Obviously, if we were using this algorithm to seek the 
as-yet unknown solution to a 3-SAT problem, we would 
not know in advance whether this problem would prove 
to be at the ‘fast’ or ‘slow’ end of the spectrum. (In¬ 
deed, we would not know whether there is a solution and 
whether it is unique; I consider this point later, but sup¬ 
pose for the moment that it is known that any solution 





































































is indeed unique.) We would then need to make a choice 
of how many clause checks to require before calling a run 
a success and measuring the qubits, as called for by the 
overarching algorithm Fig. 4. I will use the symbol A^fuii 
for this user-selected number. A simple strategy would 
be just to choose iVfuu to large enough that the bulk of 
the distribution of A^hiRes, say 99.9% of it, is lower. For¬ 
tunately there is relatively little cost in setting a high 
value for A^fuii, because once a high fidelity approxima¬ 
tion to |0) is reached then of course additional clause 
checks are unlikely to fail (but see the effect of imper¬ 
fectly applied operations, discussed later). As for the 
possibility that the 3-SAT problem in question is an ex¬ 
ceptionally ‘slow’ one, beyond even our supposed 99.9% 
threshold, assessing the impact of this lies beyond the 



Nf^iFid/ the number of clause check 
cycles until fidelity target reached 

FIG. 5. For a given number of booleans, a randomly selected 
3-SAT instance (with a unique solution) may correspond to a 
relatively quickly evolving quantum system or a rather slower 
one - i.e. there is variation in the rate at which the quantum 
state evolves to a target state |0) according to the specific 
nature of the problem . I have not studied the dependence 
in detail but, for example, having multiple solutions near the 
unique correct solution tends to slow it down. As the prob¬ 
lem size (number of booleans) increases then this distribution 
changes - but over the range I can inspect the distribution 
spreads only modestly. 

The right-side tails of these distributions appear to be ‘heav¬ 
ier’ than exp(—fcA'^), and are reasonably fitted by exp{—kN). 


range of the simulations and analysis I present here - it 
may be that the distribution of A^hiRes has a hard upper 
limit so that this problem does not arise; conversely it 
may be the case that there is a signature in the statistics 
of measured qubits which can indicate a solution exists 
but that Affuii should be increased in order to find it. 

The ‘sculpting’ algorithm shown in Fig. 4 calls for us to 
repeatedly create an approximation to |0) and measure 
it. Each successful create-and-measure process will yield 
information that is correlated to the 3-SAT solution (as 
discussed later). But what is the expected number of 
clause checks Ctotai needed for one successful full run? 
We can write 

Ctotai = Rrps fllcPR 

where i?Rps is the expected number of runs that must be 
attempted in order to achieve one successful run (“RPS” 
for Runs Per Successful run), and Ccpr is the expected 
number of clause check operations per run (“CPR” for 
Clauses Per Run). Now the former is simply Rpps = 
(Psucceed)”^ where Succeed IS the probability of a run suc¬ 
ceeding, and as noted earlier this is well approximated by 
Tsucceed ~ (cos |)^" for H qubits. The latter term can be 
written as 

f^CPR “ A^full Psucceed T ^ ^ t Tfail(^) 
i=l 

where Nf^n was introduced above as the user-selected 
number of clause checks in a successful complete run, and 
Ptniiii) is the probability that a run aborts specifically 
during the clause check. Therefore 

Ctotai = A^fuii + F f sec^ - j where F =^i pfaii(f). 

(4) 

As noted in Fig. 5 and its caption, over the limited 
range for which I have performed numerics, the suitable 
choice for A^tuu appears to scale only linearly with n, the 
number of booleans.^ Taking this to be the case, then 
the scaling for F must be at-most linear, since it is the 
expected number of clause checks in a run, given that 
the run fails i.e. it does not achieve Nfuu clause checks. 
Figure 6 shows how F behaves over a range of problem 
sizes. One sees a linear dependence becoming sub-linear 
for large n, and even declining as n becomes very large. 

If the trends appearing in the numerics of Fig. 5 and 
Fig. 6 prove to continue to large n, in particular if the re¬ 
quired A^fuii continues to scale only linearly with n, then 


® In fact Alfull must scale at least linearly, since we cannot create 
the solution-encoding state |0) without having evaluated every 
clause at least once, and for non-trivial 3-SATs the number of 
clauses scales linearly with n. 
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this algorithm constitutes a means of creating a high fi¬ 
delity approximation to state |0) in a time scaling as 
(sec^d)) where we can choose 9 in the range 0 < 0 < 
But what use is a single instance of \6) in solving the 3- 
SAT? Recall that if 6* = | then the state \6) can simply 
be measured to reveal the (unique) solution; but of course 
then sec^(|) = 2 and we have merely matched the scal¬ 
ing of the time cost for basic random classical search. In 
general if we measure the qubits in a perfectly-prepared 
state |0) in the z—basis then each qubit (independently 






Below: 

cumulative success probability 
versus number of clause checks 
for each problem size n, showing 
that for larger n the failure events 
are moreagressively'front loaded'. 




FIG. 6. Values of F in Eqn. 4 determined from numerical 
experiments. Two curves for three solution cases are shown. 


of the others) gives the ‘correct’ truth value with prob¬ 
ability |(1 -I- sin0). Thus by repeatedly generating, and 
measuring out, the state |0) we can rapidly reach a high 
degree of certainty as to the 3-SAT solution. An anal¬ 
ysis is presented in Appendix 4; in short, the expected 
number of repetitions R needed to infer the full 3-SAT 
solution grows only as a logarithmically with n. 

Now in practice we cannot create a perfect |0) in finite 
time for 9 < ^, however as shown in Fig. 5 we can reach 
a high fidelity for the entire state, and then each qubit’s 
probability of correlating with the correct truth value will 
indeed approach the ideal of |(1 -I- sin0). 

A hybrid approach 

The previous two sections have described two ap¬ 
proaches to measurement-driven 3-SAT solving: an 
‘adiabatic-like’ approach where the key parameter 9 is 
slowly evolved, and a ‘sculpting’ approach where 9 — 9o 
is fixed and consequently the system’s state changes dra¬ 
matically at first and then is more slowly honed. A draw¬ 
back of the latter is that we need to repeatedly create- 
and-measure \ip) in order to glean enough information to 
solve the 3-SAT. 

A third possibility is to create the state \6) for some 
^0 < f via the sculpting approach, and then continue 
with a slow evolution to 0 = f by the adiabatic-like 
approach (see Fig. 7). Then we gain the entire 3-SAT 




FIG. 7. A hybrid algorithm. Instead of the repeating prepare- 
and-measure process in the sculpting algorithm (upper panel 
here) we need only ‘sculpt’ |0) once, using some intermediate 
value of 0 = 00, and we proceed to evolve slowly to 0 = 
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solution by measuring the register once, at the end of 
the first successful run of the algorithm. 

Whether or not this approach is workable depends on 
the cost of the second phase: inevitably there will be 
some probability of failure during this phase, and (as we 
saw from the analysis of the original adiabatic-like ap¬ 
proach) the cost of this failure is not necessarily reduced 
when the evolution is slowed. However, at least for some 
values of 0o the hybrid approach does indeed succeed, as 
shown in Fig. 8 and Fig. 9. 

In Figure 9 three measurement-driven approaches com¬ 
pared. Table IV summaries the performance, with the 
caveat that there are further optimisations possible for 
the ‘hybrid’ approach (I made the arbitrary choice to 
split the evolution ‘50/50’ between the static 6 phase and 
the ramping phase, and to use a linear ramp). 

For each approach, the table lists under ‘Cycles’ the 
number of complete clause check cycles that were found 
to be optimal in solving a particular 24 boolean 3-SAT 
problem (the same problem used in Fig. 2). This is 



0 100 200 300 400 500 



FIG. 8. Hybrid algorithm with 00 = f- By moving sufficiently 
slowly, the failure probability in the second phase is made very 
small. 


the user-selected number that determines how slowly 0 
evolves to However most attempts will fail; conse¬ 
quently we expect to have to try multiple times and this 
is listed under ‘Tries’. But the important figure is the 
one in the rightmost column: how many clause checks do 
we expect to have to perform in total before we eventu¬ 
ally have a successful run and therefore find the 3-SAT 
solution. 




FIG. 9. A comparison of the three variations on the 
measurement-driven algorithm used to solve a particular 24 
boolean USA 3-SAT problem which was analysed earlier 
in Figs. 2 & 3. The green curve here corresponds to the 
green curves in those earlier figures. The parameters and 
performance of these three approaches are summarised in 
Table IV. The best performance is obtained using the hybrid 
approach, i.e. the ‘stepped 6’ which starts from 6 — (0.56)^ 
and holds this fixed for 37 full clause cycles before linearly 
ramping to 9 = ^ over the remaining cycles. 
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TABLE IV. Comparison of approaches 


Approach 

Used 

Optimal 

Cycles 

Tries 

Needed 

Total 

Time 

Grover’s Algorithm 

N/A 

1.18"- 

288,252 

Adiabatic-like: Linear ramp 

92 

83 

386,462 

Adiabatic-like: Sqr-root ramp 

73 

99 

190,885 

Hybrid with ‘stepped’ 9 ramp*’ 

75=37-h38 

143 

38,031 


Tries is 1.18 rather than 1, see Appendix on Grover’s 
Using 9o = (0.56)7r/2 which is optimal for this problem. 


approach would be as follows: In a multi-solution 3-SAT 
problem whenever a given boolean varies between the 
solutions (i.e. it is TRUE for some solution(s) but not 
all solutions) then the corresponding qubits will show a 
lesser asymmetry (frequency of O’s versus I’s) as com¬ 
pared to the expected bias (1 -I- sind)/2. Noticing this, 
one would reduce the SAT problem by setting the corre¬ 
sponding boolean to TRUE, say, in order to narrow the 
search for a satisfying solution. This would be repeated 
until the qubits in the reduced problem show the proper 
correlation to unique truth vales. 


VI. MULTIPLE 3-SAT SOLUTIONS 

To this point, I have considered only 3-SAT problems 
that have a unique satisfying assignment (USA). The ex¬ 
tension to cases with multiple solutions is straightfor¬ 
ward. For the USA case, for a given 9 there was a single 
state that passed all the stabiliser checks with certainty, 
i.e. Eqn. (2). Now instead there is a subspace spanned 
by a set of such states, one for each solution of the 3-SAT. 
These states are not orthogonal (except when d = ^), but 
they are mutually independent. 

The goal of the algorithm will now be to drive the qubit 
register into this subspace. There will be no guarantee 
as to whether the final state contains equal amplitudes of 
the different solutions, and indeed in general it will not; 
however provided that we can learn one of the solutions 
then we can modify the 3-SAT problem to exclude this 
solution (for example by adding a clause which the found 
solution does not satisfy) and thus search for others.^ 

One can obtain a fidelity metric to measure how well 
our qubit register is driven towards the satisfying sub¬ 
space (see Appendix 2), and thus we can ask how long 
it takes to reach the same ^ 0.999 pre-qubit fidelity cri¬ 
terion that we used earlier for the USA problems. Of 
course, now only the qubits that correspond to booleans 
whose value is the same over all solutions will approach a 
specific state, and as we might expect the algorithm runs 
a little faster. A couple of such plots are presented in 6, 
for comparison with unique solution cases. 

Obtaining one full solution is straightforward if we are 
taking the simple adiabatic-like approach (Fig. 1) or the 
hybrid ‘two phase’ approach (Fig 8) which both conclude 
with 9 = ^, since then merely measuring all qubits in the 
z-basis will reveal a solution. 

If instead we are taking the sculpting route (Fig. 4) 
of preparing \6) for some 0 > 0 and then measuring for 
partial information, before repeating for additional in¬ 
formation, then things are a little more complex. One 


^ A difficult case would be a 3-SAT with two solutions differing in 
only one boolean, so that adding a clause to exclude the found so¬ 
lution would likely exclude the alternative solution also. In that 
case Order(n/3) further experiments would suffice to find the 
second solution, by trying a series of different additional clauses. 


VII. ROBUSTNESS VERSUS ERRORS 

There is reason to hope that the approach described 
here can have a degree of inherent error tolerance, where 





FIG. 10. The effect of adding noise to the rotations during 
the running of the ‘sculpting’ algorithm (6 = (0.5)^) on a 
particular 20 boolean problem. Note that noise can boost 
the fidelity (the green line goes above the blue in places). 
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by ‘error’ I mean imperfections in the preparation, ro¬ 
tation, and measurement of qubits, or general decoher¬ 
ence to the qubits from the environment. Because the 
approach works by continually projecting the system to¬ 
ward a target state |0), certain small deviations to the 
state may be corrected. Meanwhile large deviations may 
cause the protocol to abort, but that is also a form of 
error handling. 

As a preliminary exploration of this possibility, I intro¬ 
duced a small random contribution to every single-qubit 
rotation in the sculpting algorithm, i.e. each rotation 
over or under rotates by a random defect that is capped 
at a certain percent. The results are shown in Fig. 10. 
One can see that indeed the algorithm still functions, 
albeit at a slightly reduced level. 

A proper exploration of error tolerance would of course 
require far more that this - general single-qubit er¬ 
rors (rather than merely over/under rotating around the 
y—axis) and moreover it would be important to consider 
correlated errors from the three-qubit gate and ancilla 
measurement. 


VIII. CONCLUDING REMARKS 

The numerical simulations cannot establish whether 
this approach would perform well in the regimes where 
classical approaches struggle; possibly analytical consid¬ 
erations may be able to rule out this promising prospect 
(assuming that further numerics do not already reveal 


some error in the data presented here). One point of 
concern arises from Fig. 5 where we see that certain 3- 
SAT problems are ‘harder’ for the quantum system than 
others - these harder cases evolve more slowly towards 
the target \6) state. I do not have a complete expla¬ 
nation for this, although I have noted that the number 
of nearly successful boolean assignments is correlated to 
how slowly the state evolves, as one might expect. 

My original motivation for exploring these algorithms 
was to construct a circuit based approach to optimisation 
that has some inherent error tolerance, it a way analogous 
to adiabatic quantum computing. I am optimistic that 
this may indeed prove to be the case, although the pre¬ 
liminary simulations I have performed so far are not so¬ 
phisticated enough to incorporate general noise sources, 
only the highly specific case of over/under rotation. But 
for that special case, the tolerance seems good. 
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Appendix 1: Grover’s algorithm 

Regarding the use of Grover’s algorithm to directly 
solve a 3-SAT of 24 booleans: we have N = 2^^ can¬ 


didate solutions given the 3-SAT has a unique solution. 
We could set Grover’s to perform m iterates such that 
sin^((2m -I- 1)/'/N) is close to unity, which would give 
us m ~ 2^°7r ~ 3,217 iterates. However because of the 
sine squared function, the last ^ 20% of the iterates give 
us a relatively poor return on the invested time. Con¬ 
sequently the optimal approach is “stop early” at 2,386 
iterations and measure; at that point the probability that 
our measurement gives us the correct solutions is ~ 84% 
but we simply restart the measurement does not give us 
our solution. The expected number of ‘runs’ needed is 
then 1.184. Then the expected total number of iterates 
to reach a solution is reduced to 2, 826, with each it¬ 
erate requiring all the clauses to be tested. There are 
round(4.267* 24) = 102 causes, so that the expected total 
number of clause checks required by Grover’s algorithm 
is 288, 252. However it is worth reiterating that this ap¬ 
proach requires this expected number of clause checks for 
any 24 boolean 3-SAT with a unique satisfying solution - 
i.e. the particular features of the problem are irrelevant. 

Appendix 2: Fidelity is monotonic 

This section considers the evolution of the qubit reg¬ 
ister’s fidelity with respect to the target state, for a 3- 
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SAT with a unique satisfying argument. However the 
extension to the multi-solution fidelity (next section) is 
st raight forward. 

As observed by van Dam [7], for certain adiabatic algo¬ 
rithms one can construct 3-SAT problems for which the 
evolution ‘gets stuck’ and never finds the actual solution. 
Generally this is a potential issue in any approach (quan¬ 
tum or classical) that seeks a solution by minimising a 
cost function. An example of a cost function is simply 
the number of unsatisfied 3-SAT clauses. An approach 
that involves updating the system’s state so as to follow 
the cost function ‘down hill’ to reach the lowest possible 
cost, may run the risk of being trapped in a local mini¬ 
mum (and therefore will require a strategy for escaping 
such minima). 

The present approach does not appear to have this 
issue: The fidelity of the qubit register with respect to 
the target state will only increase with each complete 
cycle of the clause checks. (This is not unexpected since 
the obvious quantum approach, i.e. Grover’s algorithm, 
also does not involve a cost function and cannot become 
trapped.) Generally we can write the state of our register 
as 

\^) = a\e) + I3\e^) 

where 6 is the desired state as 6^ is some state orthog¬ 
onal to the desired state. The fidelity is |ap and pro¬ 
viding that this is non-zero then when the next clause 
check operation is performed there is a finite probability 
of ‘passing’ that check; then the post-check state is 

\^) = J^{a\e) + I3'\e^)) 

where /3' < /?, normalisation constant Af > 1. 

The inequalities are not strict because we have not ex¬ 
cluded the possibility that the three qubits involved were 
certain to pass the clause check - they could have been in 
a superposition of several states, all of which satisfy the 
clause in question but only one of which satisfies the 3- 
SAT. However a complete cycle of successful clause check 
operations must indeed reduce /3, unless it was already 
zero, because (as discussed earlier) |0) is precisely the 
unique state that passes all clause checks with certainty. 

Of course this leaves open the question of whether one 
can adversarially construct a 3-SAT problem that has an 
extremely slow evolution toward the target state |0). 

Appendix 3: 

Fidelity given mnltiple solntions 

If we have a 3-SAT with m > 1 solutions, then applying 
any of the algorithms discussed here will drive the qubit 
register toward a subspace spanned by the states 

\0j) = ^ y {Lij9) |-|-)i for j = 1... m (5) 

i 

where the various sets of values {L}j correspond to the 
different 3-SAT solutions in the manner specified after 


Eqn. (2). Note that the value of 9 is of course the same 
for all states; only the L values vary. 

We can write a general state of our n qubit register as 

I'f') = l^i) + |iA) 

where |V||) lies in the solution subspace, and |VD is 
orthogonal to it. The probability |(V|||V||)p is a mea¬ 
sure of the fidelity with which I'k) approximates an ideal 
solution-state, i.e. any state completely within the sub¬ 
space. We may wish to record and plot this quantity 
during a simulation of an algorithm. 

By definition 

(6>j|Vl)=0 and |V||) = ^ c,|6/,) 

i 

remembering however that the \6j) are not mutually or¬ 
thogonal. Then 

(0j|'k)=^c,(0j|0,) or v = Mc 

i 

where vector v has elements (6j|'I'), vector c has ele¬ 
ments Ci, and M is a matrix of elements mij = {6j\6i). 
This matrix is easy to compute since the defition of the 
solution states Eqn. (5) implies that projection {6j\6i) 
is simply (cos 9Y where p is the number of mismatch¬ 
ing booleans (i.e. the number of instances where a given 
boolean is TRUE for solution i but FALSE for solution j, 
or vice versa). 

It is useful to compute the inverse matrix be¬ 

fore beginning the simulation of an algorithm. Then at 
any time during the simulated execution, when the qubit 
register is in some state j^*), we can determine the Ci 
by finding the elements of vector v, i.e. the projections 
(0i|'I'), and simply evaluating 

c = M-^ V. 

Thus we can construct |y||) and moreover we can com¬ 
pute and record the projection on the desired subspace. 

Appendix 4: 

Information from measuring a 0 < | state 

Suppose that a high fidelity approximation to the state 
|0) with 0 < f has been created by e.g. the sculpting 
algorithm. This state is simply a product of individ¬ 
ual qubit states F(±0)|-|-) where the sign is positive if 
the corresponding boolean is FALSE in the unique sat¬ 
isfying assignment (USA), and negative otherwise. If 
a given qubit is measured in the z-basis then the mea¬ 
surement will give the correct truth value (according to 
rule |0) —FALSE and |1) —>■ TRUE) with a probability 
p = {1 + sin6*)/2. In the following I will call a measure¬ 
ment ‘incorrect’ if it implies the wrong boolean value; 
note that I do not mean to suggest some failure in the 
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measurement apparatus (which I am taking to be per¬ 
fect). Measurements are independent of one another in 
this sense: the question of whether a given qubit gives 
the correct value is not correlated to the correctness of 
any another qubit measurement. (Note that here we are 
considering the exact state |0), whereas of course the real 
state created using e.g. the sculpting protocols will not 
be exact and may include highly correlated elements - 
but for a sufficiently high hdelity the remarks here will 
be valid.) 

Having measured all qubits and noted the results, we 
can repeat the process of creating \6) and again measure 
it. This new set of measurements will be independent 
of the former set. So the probability of a measurement 
being correct is simply p for every qubit and on every 
run. 

The interesting question is, how many runs will it take 
until we can correctly guess the boolean truth values? 
We know the value of p but we do not know whether a 
given qubit is biased towards measurement outcome |0) 
(implying the boolean is FALSE) or outcome |1) (imply¬ 
ing TRUE). Assuming we always measure in the z-basis, 
then this scenario with an n booleans is equivalent to the 
following classical problem: we are given n numbered 
coins, each with the same known bias p > ^, but we do 
not know whether a given coin favours ‘heads’ or ‘tails’. 
We are allowed to throw all coins at once and note the 
outcomes; how many such throws do we need before we 
can guess the biases correctly? 

There are several ways to make this question precise 
and then derive an answer. I prefer the following: What 
is the expected number of throws, i?, required until our 
best guess of the set of boolean values is probability er¬ 
ror free - i.e. the guess is fully correct with probability 
greater than a half. Now after R throws we make our best 
guess about the bias of a given coin, i.e. head-favouring 
or tail-favouring, simply by noting whether it produced 
more heads or more tails (we can assume R is odd to 
avoid even-splits where we don’t know how to guess). 
Let’s say a given coin is head-favouring. Then the prob¬ 
ability that we wrongly guess ‘tail-favouring!’ is given by 
the sum of the probabilities that the number of heads, i, 
is less than half of R. 


will probably not occur over the entire set of n coins, we 
require that Pwrong assuming that n is large, i.e. 


i=M 

E 

i=0 


^ I p\l-p)^-" < -. 

i n 


If the product R{1 — p) is greater than about 5, the 
sum is well approximated by an integral. Making that 
assumption, we have: 




^ {x-Rpf ^ 


dx < — 
n 


where the variance cr^ = Rp{l ~p)- Rescaling x yields 

n-G 


ip - 

We recognise the complimentary error function erfc. 


exp(-x^)d.<l where 

J-oo n ^2p{l-p) 


letfc(G) < i 

and recall that to lowest order (for large x) 


erfc(a;) 


exp(—x^) 


so that the criterion becomes 
Gexp(G^) > 

ZyJ TT 

Now returning from coins to qubits and noting the as- 
sumptions made, we conclude: When attempting to solve 
USA 3-SAT problems of n booleans using the sculpt¬ 
ing algorithm with R repetitions of the “prepare |0) and 
measure” cycle, R need only inerease logarithmically with 
problem size n. 

Moreover we recall that bias p = (1 -I- sin0)/2 so we 
can rewrite 


sin9'/R 

-I- sin6*)(l — sind) 



tan0. 


Pwrong = E ^ ^ ^ \ 

where M = {R — \)/2. In order that such a mistake 


We might also consider a 3-SAT of a hxed size and 
ask how the selected value of 9 influences the expected 
number of repetitions R needed to guess all booleans cor¬ 
rectly. This can be assessed simply by requiring G to stay 
constant, so Roe cot^ 9 and thus R ~ 9~‘^ for small 9. 



